Spin-Peierls Transition in CuGeOs: Critical, Tricritical or Mean Field? 
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The spin-Peierls phase transition in CuGeOs has been extensively studied utilizing a variety of 
experimental techniques. Interpretations of the phase transition behavior vary from tricritical to 
mean field to Ising critical to XY critical. We show that the behavior in the vicinity of the phase 
transition of each of the order parameter, the magnetic energy gap and the heat capacity can 
be quantitatively fitted with few adjustable parameters with a mean field model incorporating a 
tricritical to mean field critical crossover in the transition region. 

PACS numbers: 75.40.Cx, 75.40.-s, 75.30.Kz, 64.60.Fr 



I. INTRODUCTION 

The spin-Peierls transition corresponds to the dimer- 
ization of a one-dimensional S = h antiferromagnetic 
chain coupled to a three dimensional elastic medium jjj- 
H . Until relatively recently, spin-Peierls transitions had 
only been observed in organic charge transfer compounds 
such as copper bisdithiolene (TTF-CuBDT) §- g. Ex- 
perimental information obtainable in such systems has 
been limited both by the size of available single crystals 
and by the sensitivity of these materials to damage by 
x-rays or electrons. Nevertheless, some important infor- 
mation on the spin-Peierls phase transition has been ob- 
tained in a number of different organic materials. Inter- 
estingly, in most, if not all cases, the data are consistent 
with a simple BCS-type mean field transition S- ||. 

Much more complete experimental work on the spin- 
Peierls transition has been made possible by the discov- 
ery that a structurally simple, inorganic chain compound 
copper germanate (CuGeOs) undergoes a spin-Peierls 
transition at a transition temperature around 14K 
The crystal structure of CuGeC>3 is orthorhombic, space 
group Pbmm, with a unit cell of dimensions a = 4.81 A, 
b = 8.47 A and c = 2.94 A at room temperature [jl0|. 
The Cu 2+ ion carries a spin S = h and forms a (CUO2) 
chain with the neighboring Cu 2+ ions along the c-axis 
direction. The successive Cu 2+ S = | spins are antifcr- 
romagnetically coupled through the superexchange inter- 
actions via the bridging oxygen atoms. Below the spin- 
Peierls transition temperature, T$p, the dimerization of 
Cu-Cu pairs along the c-axis direction, accompanied by 
shifts of the bridging oxygen atoms in the ab plane, gives 
rise to superlattice reflections at the (|, k, I) ( h, I: odd 
and k: integer) reciprocal-lattice positions |ll|]. These 
have been observed in electron diffraction (l2) , x-ray [jl3| , 
and elastic neutron scattering [Tl| ] experiments. Using 
coarse resolution x-ray diffraction techniques, Pouget et 
al. |D| have measured the pretransitional thermal lattice 
fluctuations whose correlation lengths diverge anisotropi- 



cally with decreasing temperature in a manner consistent 
with mean field theory. These same fluctuations have 
been studied at high resolution using synchrotron x-ray 
diffraction techniques by Harris et al. 14 1. These lat- 
ter authors observe within about IK of T$p large length 
scale fluctuations with characteristic length scales about 
an order of magnitude longer than those characterizing 
the bulk critical fluctuations. 

In spite of this large amount of work, it is still not 
agreed whether the observed transition behavior reflects 
mean field or critical behavior. Extant models include: 
a) tricritical to 3D Ising crossover behavior |l4||nj; b) 
mean field behavior |^5|; c) 3D XY with corrections to 
scaling |l6| , and, most exotically, d) a 2D XY to 3D XY 
crossover as T$p is approached ]l7| ]. Harris et al. |l4| ] 
first argued that because of the one-component nature 
of the dimerization order parameter for a spin-Peierls 
phase transition, asymptotically the transition must be 
in the 3D Ising universality class. They argued fur- 
ther, that because of the coupling to the elastic strains, 
the precritical behavior should be tricritical-like. Sim- 
ilar conclusions, albeit based on different physical rea- 
soning, were arrived at later by Werner and Gros [fT5"| . 
Proponents of 3D XY behavior typically argue that the 
copper and oxygen displacements are independent thence 
yielding a two-component order parameter Jl8| . Implic- 
itly, Harris et al. Jl4| assume that all of the atomic dis- 
placements accompanying the spin-Peierls transition are 
linearly coupled thence reducing the system to a one- 
component order parameter. The 3D critical behavior 
models seem to be supported by measurements of the 
order parameter |]U4] , |l6| , p^7| which for reduced tempera- 
tures - 2 x 1CT 3 < 1 - T/T S p < 5 x 1CT 2 exhibits 
power law behavior (1 - T/Tsp) 13 with /3 = 0.33 ± 0.02, 
in good agreement with both 3D Ising and XY values of 
(i = 0.325 and 0.345 respectively ]19j . The heat capacity 
data are equally well described by a 3D critical behav- 
ior model (Ising or XY) and by a mean field model with 
Gaussian fluctuations lEOfl. 
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In this paper we present an alternative model for 
CuGeC>3, namely a Landau- Ginzburg model incorporat- 
ing a tricritical to mean field crossover. As we shall show, 
this model describes all available data very well with few 
adjustable parameters. The format of this paper is as fol- 
lows. In Section II we introduce the model including its 
genesis in studies of critical phenomena in thermotropic 
liquid crystals systems. Section III presents an analysis 
of the available data for CuGe03 using this model. Fi- 
nally, in Section IV we give a summary, our conclusions 
and suggestions for future experiments. 



where Tcr is the crossover temperature from tricritical 
to mean field behavior. We note that in the above expres- 
sions the exponents are fixed and only the amplitudes and 
the two temperatures, Tsp and Tcr, are variable. A log- 
log plot of Eq. (6) reveals that for the order parameter <j> 
the effective exponent (3 crosses over gradually from j to 
\ as T varies from less than to greater than Tcr- In the 
smectic A - smectic C case the measurements span Tcr 
and accordingly intermediate exponents, ~ 0.36, are 
found even though the actual transition is mean-field- 
like for temperatures in the immediate vicinity of Tac 



II. THE MODEL 

The conundrum described above for CuGe03 is remi- 
niscent of a similar divergence of views which occurred in 
the interpretation of experiments on smectic A - smectic 
C phase transitions in thermotropic liquid crystal sys- 
tems |2lj- |23| . In particular, in that case, measure- 
ments of the tilt order parameter |2l],|23j typically reveal 
power law behavior (f> ~ (1 — T /Tac) ^over the temper- 
ature range 5 x 1CT 5 < (1 - T /Tac) < 5 x 1CT 3 with 
(3 = 0.36 ± 0.02. However, this divergence of views was 
resolved by Huang and Viner |22] ] and Birgeneau et al. 
p3"t who showed that all of the data including the heat 
capacity, order parameter, and tilt susceptibility, were 
consistent with the predictions of a simple Landau model 
with an anomalously large 6th order term. Clearly, it is 
of interest to carry out a similar analysis for the available 
data for the spin-Peierls transition in CuGeOa. 

For the Landau-Ginzburg model the free energy is 
given by 



F = ar<j) 2 +b<p i + C(f) 6 .... + 
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where r = T/T c - 1. 

With To = b 2 /ac, standard calculations yield for the 
order parameter, <j), the specific heat, C, the susceptibil- 
ity, Xi an d the correlation length, 



C = 



0=(V3c) 1/2 [(l-3r/r o ) 1 / 2 -l] 1 / 2 
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with similar expressions for r < for \ an d £■ Eq. (2) 
and (3) are conveniently rewritten in the form 
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III. ANALYSIS 

We now apply this tricritical-mean field crossover 
model to CuGeOs. The first test is T$p itself or, more 
precisely, the ratio of the spin gap, A, to Tsp- In the 
mean field theory of Pytte |24j| , the spin-Peierls tran- 
sition is BCS-like so that in the weak coupling limit 
2 A /T S p = 3.5. In the charge transfer salts TTF - 
CuBDT TTF - AuBDT @, MEM - (TCNQ) 2 

§, and SBTTF - TCNQC1 2 § this ratio is found to 
be 3.5, 3.7, 3.1 and < 3.5 respectively, in good agree- 
ment with the BCS value. Critical fluctuations, either 
Ising or XY in character, would act to increase this ra- 
tio. For CuGeC-3, A = 24.5K and T S p ~ 14K implying 
2A/Tgp — 3.5, consistent with a BCS mean field the- 
ory description P,|25[|. At the minimum, this value for 
2 A /Tsp argues against any quantitatively important ef- 
fect of critical fluctuations on Tsp in CuGe03. 

The behavior of the order parameter in CuGeOa is 
of particular importance since this observable appears 
to provide the strongest evidence for true critical rather 
than mean field or tricritical behavior. A number of 
groups have reported measurements of the temperature 
dependence of the order parameter in CuGe03 Jl4|,[l6|,|l7j . 
The measured phase transition temperature T$p varies 
between 13.3K and 14.6K in different samples. Never- 
theless, near-universal behavior is observed for the order 
parameter provided that it is plotted as a function of 
the reduced temperature T/Tsp- As noted above, fits 
of the order parameter <p(T/T SP ) for 1 - T/T SP < 0.05 
to a single power law <jf> ~ (1 — T/Tsp)@ all yield val- 
ues of (3 = 0.33 ± 0.02. As discussed by Gaulin and 
co-workers jl6f| , inclusion of a correction-to-scaling mul- 
tiplicative factor (1 + -BItI" 5 ) in the expression for <p both 
improves the goodness of fit and, not surprisingly, ex- 
tends the range of validity of the fit. 
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FIG. 1. (3/2, 1, 3/2) superlattice peak intensity mea- 
sured with synchrotron x-ray diffraction techniques. The 
peak intensity is proportional to the order parameter squared, 
<h 2 . The dashed line is the result of a fit of the data for 
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i — x — < 0.04 to a single power law <f> 2 ~ cAq(1 ~ t^T" > 
with p = S 0.314 ± 0.01. The solid line is the result of alft to 
Eq. (6) with t C r = 0.006. 



We show in Fig. 1 our own measurements of the 
order parameter squared in a sample of CuGeC>3 with 
Tsp — 14. 6K. These data are consistent with those mea- 
sured by both ourselves and other groups in a variety of 
samples |l4] , |l6| , p7| . Fits to a single power law for r < 0.04 
yield /3 = 0.3 14 ±0.01. However, as noted by Harris et 
al. JIJ] and as may be seen in Fig. 1, the data fall sig- 
nificantly below the power law curve for r > 0.04. We 
show, in addition, in Fig. 1 the results of a fit to the 
tricritical to mean field crossover form, Eq. (6). This 
fit has only three adjustable parameters, c/>q, Tcr and 
Tsp. This is the same number of parameters as those 
in the single power law fits discussed above and two less 
than the number of adjustable parameters in fits to a 
power law with corrections-to-scaling with both B and 8 
varied. (We note that Lumsden et al. |l(| fix 6 = 1/2 
whereas Lorenzo et al. allow 6 to vary; the latter group 
find an optimum fit for 5 ~ 1). It is evident that Eq. 
(6) describes the order parameter data extremely well 
over the complete range of temperatures. The fit yields 
T CR = 1 - Tcr/Tsp = 0.006 ± 0.001 implying that the 
crossover from tricritical to mean field behavior occurs 
at a quite small reduced temperature. 

We now discuss the energy gap A. Using a simple 
scaling ansatz, Cross and Fisher |3| argue that A ~ c/> 2 / 3 . 
We show in Fig. 2 the data of Lorenzo et al. |l7j for the 
magnetic energy gap for T < T$p in a sample of CuGeC>3 
with Tsp = 14. 4K. In part because of the apparent jump 
of A(T) at Tsp, Lorenzo et al. Jl?]] interpret these data as 
indicating a 2D XY Kosterlitz-Thouless transition p6| . 
In fact, these data are readily explained using the model 



of Cross and Fisher |J together with the tricritical-mean 
field crossover form for <f>, Eq. (6). In this case we hold 
T SP fixed at T SP = 14. 4K and set t cr = 0.006 as deter- 
mined above so that there is only one adjustable param- 
eter, the overall amplitude A(0). The result so-obtained 
is shown in Fig. 2. It is evident that the tricritical-mean 
field model with A(T) ~ c/> 2 / 3 describes the measured gap 
energy A(T) extremely well over a wide range of temper- 
atures with only one adjustable parameter. Indeed, this 
is by far the best test to-date of the Cross-Fisher model. 
We should note that this model cannot explain the in- 
ferred pseudogap above Tsp [jl7|. However, the "pseu- 
dogap" is deduced using a heuristic line-shape analysis 
which lacks a firm theoretical basis. 
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FIG. 2. Magnetic energy gap in CuGeOa. These data are 
from Ref. 17. The solid line is the result of a fit to the form 
A(T) = A(0)c/> 2/3 where cj> is given by Eq. (6) with tcr held 
fixed at 0.006. 

The specific heat in CuGeOa has proven to be the 
most difficult thermodynamic quantity to interpret un- 
ambiguously [^0|. This is, in part, because of the extreme 
sensitivity of the specific heat near Ts p to sample inho- 
mogeneities and, in part, because of the inevitable large 
number of adjustable parameters required to describe the 
critical specific heat in any physically relevant model. 
Fig. 3 shows high resolution magnetic specific heat (Cm) 
data for a sample of CuGe03 with Tsp = 14.24K from 
Lasjaunias and coworkers (20). Hegman et al. have 
carried out an extensive analysis of these data using both 
a mean field "BCS plus Gaussian fluctuation" model and 
a critical behavior model. They find that both mod- 
els describe Cm quite well in the immediate vicinity of 
Tsp, albeit at the cost of a rather large number of ad- 
justable parameters. The critical behavior model fits give 
a value for the specific heat exponent, a, near 0. On 
the other hand, the Gaussian fluctuation analysis implies 
that the true critical behavior is confined to the region 
|t| < 0.0006. 
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FIG. 3. Magnetic specific heat in CuGeC>3. These data are 
from Ref. 20. The solid line is the result of a fit to Eq. (8) 
with tcr held fixed at 0.006. 

Given the uncertainties connected with the fits de- 
scribed above, the best one can hope for is to determine 
whether or not the tricritical-mean field crossover model 
is consistent with the experimental results for Cm shown 
in Fig. 3. First, it is evident that Eq. (7) will be inade- 
quate since one must, at the minimum, include Gaussian 
fluctuations above Tgp. We therefore include the fluctua- 
tions above T$p in the simplest way possible by replacing 
Eq. (7) by 



Cm 



C M T (l + 3 T J p 5f£ R 



C~ M T ( 1 



Tsp—Tcr 



-1/2 
-1/2 



■ 7T r > 

B_ t < 

(8) 



where 7T is the regular linear term for a ID Heisenberg 
antiferromagnet and B_ is the background term below 
Tsp- The background _B_ should, in general, be tem- 
perature dependent; however, given the narrow range of 
temperatures we consider, a constant background is ade- 
quate. Eq. (8) is closely similar to the BCS plus Gaussian 
fluctuation model considered by Hegman et al. f20(| since 
the Gaussian fluctuations give rise to a |t| -1 / 2 contribu- 
tion to Cm both above and below T$p. The solid lines 
in Fig. 3 correspond to fits to Eq. (8) with tcr fixed 
at 0.006 and C M , C M ,j, B_ and T$p varied. Clearly 
Eq. (8) describes Cm quite well; indeed the fit appears 
to be better than those for either of the models tested 
by Hegman et al. pc[ . The fit shown in Fig. 3 gives 
C M /C M = 1.1 ± 0.13; this ratio is expected to be non- 
universal so it cannot be simply interpreted. We con- 
clude, therefore, that the tricritical-mean field crossover 
model describes Cm well although not uniquely so. 

Finally, we discuss the correlation length and the stag- 
gered susceptibility. Pouget et al. |13| have found that 
the correlation length over a wide temperature range fol- 
lows the behavior £ ~ (T/Tsp — l) -1 / 2 , consistent with 



mean field theory; however, the number of data points in 
their experiment near Tg p is sufficiently small that their 
results do not meaningfully differentiate between various 
theoretical models. Harris et al. jl4| have reported a high 
resolution synchrotron x-ray study of the critical fluctu- 
ations above T$p in CuGeC>3. They find pretransitional 
lattice fluctuations within IK above T$p whose length 
scale is about an order of magnitude longer than those 
characterizing the bulk thermal fluctuations. The line- 
shape of the large length scale fluctuations is consistent 
with a Lorentzian-squared form. The measured critical 
exponents are v = 0.56±0.09 and 7 = 2.0±0.3 where 7 is 
the exponent characterizing the divergence of the discon- 
nected staggered susceptibility ]27j. The mean field pre- 
dictions for these exponents are v = 1/2 and 7 = 27 = 2 
whereas for 3D Ising (XY) critical behavior one expects 
v = 0.63 (0.67) and 7 = 2.5(2.64). Thus the Harris et 
al. |lj] data favor the tricritical-mean field model but 
3D Ising or XY critical models are not excluded. Precise 
measurements of the bulk staggered susceptibility using 
neutrons should yield accurate values for v and 7 and 
this, in turn, would definitively choose between the mod- 
els. 



IV. DISCUSSION 

In summary, each of the order parameter, magnetic 
energy gap, specific heat, correlation length and discon- 
nected staggered susceptibility are well-described by a 
simple Landau-Ginzburg model exhibiting a tricritical- 
mean field crossover near Tsp. Further, the ratio of 
the energy gap to T$p is consistent with the value for a 
BCS mean-field transition. We conclude, therefore, that 
CuGeOa, in common with the organic change transfer 
salts, exhibits a mean field spin-Peierls transition for re- 
duced temperatures |r| > 0.001. 

The principal remaining issue is the microscopic ori- 
gin of the tricritical behavior. Harris et al. Jl4|] argue 
that this is caused by a diminution in the effective fourth 
order term in Eq. (1), 6^ 4 , because of coupling to the 
macroscopic strain. It also seems possible that compet- 
ing nearest and next-nearest neighbor exchange interac- 
tions along the chain could generate the tricritical insta- 
bility |l^ , pq| . Specifically, Castilla et al. |^8| argue that 
the ratio of the next nearest neighbor to nearest neigh- 
bor exchange interaction along the chain is close to the 
critical value for spontaneous formation of a magnetic 
gap independent of coupling to the lattice. Heuristically, 
it seems that this could generate tricritical behavior in 
the phase diagram. Another possible source of tricriti- 
cal behavior is competition between the Neel state and 
the spin-Peierls state, that is, competition between the 
coupling of the S = 1/2 chain to the lattice and the in- 
terchain exchange coupling. Clearly, a multidimensional 
theoretical analysis of the spin-Peierls phase diagram in- 
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eluding magnetostriction, competing intrachain exchange 
interactions together with the interchain magnetic and 
elastic coupling is required. 

Of course, the mean field behavior itself in all of these 
spin-Peierls systems is not yet well understood. In TTF- 
CuBDT there is evidence for a soft phonon at very high 
temperatures and Cross and Fisher |j) speculate that 
the precursive soft mode accounts for the large length 
scale underlying the mean field behavior. In CuGe03, 
no soft phonon at all has yet been seen. Thus, the micro- 
scopic origin of the large length scale in CuGe03 remains 
to be elucidated. 

Finally, it would be very interesting to see if the puta- 
tive nearby tricritical point could be accessed by changing 
some variable such as pressure, uniaxial stress or doping. 
Masuda et al. ]2jJ have shown that replacement of Cu by 
Mg both depresses T$p and appears to drive the spin- 
Peierls transition first order. The concomitant tricriti- 
cal point could well account for the observed tricritical- 
mean field crossover in pure CuGeC"3. We note, however, 
that the actual physics of magnetic dilution in CuGeC>3 
is quite complex since dilution introduces frustration of 
the interchain elastic interaction |50|. Replacement of 
Cu 2+ by Cd 2+ (Ref. @) or Ge 4+ by Ga 4 + (ref. §§) 
both lead to mean field behavior over quite wide tem- 
perature ranges; that is, doping with these ions moves 
CuGe03 away from the tricritical point into the pure 
mean field regime. Again, further research, both experi- 
mental and theoretical, is required to elucidate these ef- 
fects more completely. 
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